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It is shown that a minimally coupled scalar field in Brans-Dicke theory 
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' yields a non-decelerated expansion for the present universe for open, flat 
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1 Introduction 

> 

The recent extensive search for a matter field which can give rise to an acceler- 
ated expansion for the universe stems from the observational data regarding the 



O 

C\) 1 luminosity-redshift relation of the type la supernovae up to about z ~ 1 |]]. This 

matter field is called a "quintessence matter" (Q-matter for short). The most 
O ! popular candidate as a Q-matter has so far been a scalar field having a potential 

which generates a sufficeint negative pressure at the present epoch 0. Amongst 
the scalar fields considered as the Q-matter, the tracker field slowly rolling down 
its potential as proposed by Zlatev, Wang and Steinhardt [[J appears to be very 
promising. Some exotic matter like the domain walls or cosmic strings also find 
themselves amongst the possible candidates M. Unfortunately most of these fields 
work only for a spatially flat (k = 0) FRW model. Very recently, Chimento et al. 
|| showed that a combination of dissipative effects like a bulk viscous stress and 
a quintessence scalar field gives an accelerated expansion for an open universe 
(k = —1) as well. This model also provides a solution for the coincidence prob- 
lem as the ratio of the density parameters corresponding to the normal matter 
and the quintessence field asymptotically approaches a constant value. Recently 
Bertolami and Martins |J obtained an accelerated expansion for the universe in a 
modified Brans-Dicke (BD) theory by introducing a potential which is a function 
of the Brans-Dicke scalar field itself. 



The present work investigates the possibility of obtaining a non-decelerating 
(q = —ad/ (a) 2 < 0) expansion for the universe in Brans-Dicke theory with the 
help of another scalar field which is minimally coupled to gravity and serves as 
the quintessence matter. Brans-Dicke theory played a major role in the attempts 
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to solve the "graceful exit" problem in the inflationary scenario prevailing in the 
early universe. It is worthwhile to explore the possibilty of an application of the 
theory towards a reasonable explanation of the late time behaviour of the uni- 
verse as well. We find that for a negative value of the Brans-Dicke parameter u, 
the theory leads to an accelerated expansion for the universe for a spatially flat 
model. The possible solutions contain the Bertolami-Martin's solution as well. 
Furthermore, the theory also leads to at least non-decelerating expansions even 
for non flat models. It should be noted that we do not modify Brans-Dicke theory 
by introducing a self interaction of the Brans-Dicke's scalar field, but rather look 
at the possibility of a non-decelerating solution with the help of a quintessence 
field within the purview of the theory itself. 



2 Field equations and the solutions 

The Brans - Dicke theory of gravity is given by the action 

where is the BD scalar field, u is the dimensionless constant BD parameter and 
L m is the Lagrangian for all other matter fields. If we assume the matter field to 
consist of a perfect fluid and a scalar field ip as the quintessence matter, the field 
equations for a Robertson- Walker spacetime are, 

{a 2 + k) _ (p m + pip) Q d0 ucp 2 

6 a* ~ a0 + 2 0*' W 



u (d 2 + k) -(pm+P^) u4> 2 dj <j> 
a a 2 2 2 a0 0' 

Here, a is the scale factor of the Robertson- Walker metric, k the spatial curvature 
index, p m and p m are the density and the pressure of the normal matter, p^ and 
p^p are those due to the quintessence field given by 

/^ = ^ 2 + vm ^ = ^ 2 -m (4) 

where V = V(ip) is the relevant potential. The wave equation for the scalar field 
■0 reads 

h • dV(ip) 
^ + 3-^ = 5 

a dip 

and the wave equation for the Brans-Dicke scalar field is 
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The matter conservation equation, 



Pm + 3 ^(Pm + Pm) = 0, (7) 

follows from the field equations. Assuming that at the present epoch the universe 
is filled with cold matter with negligible pressure, we put p m = 0, and this 
equation integrates to 

Pm = Pi/a 3 , (8) 

where p\ is an integration constant. Since our principal interest is to find an 
accelerating power law solution for the scale factor a, we shall assume that both 
a and <p are power functions of the cosmic time t in the form 

a = ai t a , = (f)it p . (9) 

We shall look at the possibilities of consistent solutions with a%, (pi, a, /3 are con- 
stants with ax, (pi being positve definite and a > 1. These constants will be re- 
lated amongst themselves and the characteristic constants of the theory through 
the field equations. 

By combining ( 0) and ( |3[) with (^) yields the expression for ip 2 as 

= ^V~2« + (2a + af3- uf3 2 - (3 2 + P)^" 2 - ^r 3a . (10) 
The potential V can be found out from the equation ( as 



V = -\ P -±t-^ + l^- 2a + -6i[{2u + 3)09 + 3a - 1)0 

+ (2a + a(3- u(3 2 - (3 2 + /5)]t /3 ~ 2 . (11) 

The wave equation for the quintessence scalar field ip ( ||), when multiplied 
by ip, looks like 

-^ = ^ + 3^ 2 , (12) 
at a 

which readily yields a first integral and hence an expression for V as 

_ k(Pi W + Aa) t0 _ 2a 
2 a? af {(3 -2a) 

- (2a + a(3- u(3 2 - (3 2 + (3) V ' 2 - (13) 



Now we demand that the right hand sides of ( [11]) and ( [T3| ) coincide. This 
leads to the required consistency relations amongst the constants. These relations 



3 



lead to quite a few possibilities of solutions for an accelerating universe (a > 1). 
These possibilities fall into two broad classes, one in which k = and the second 
where k is different from zero. 



Case 1 : k = 

In this case, the consistency condition is 

6(4a 2 - 2a + 2a 2 (3 - a(3) = tu(3(p 2 + 6ap + 12a - 4). (14) 

It is easily seen that the condition ( [L4]) is automatically satisfied if (3 = —2, 
and so equation ( |lOj) reduces to 



2(2u + 3)0! 



3a 



This indicates that u < —3/2 as ip 2 cannot be negative. For a 
2\2uj + 3 1 0i > 2p. In this case equation ( |T5|) integrates to 

A 

* = ± T 



(15) 

4/3, one has 
(16) 



where A 2 = —2{2uj + 3)0i — ^ > 0, and a simple form of V(ip) can be obtained, 



V = Viip 4 



(17) 



Vi being a constant, related to the other constants of integrations like a\, pi etc. 
through the field equations. With a = 4/3, deceleration parameter q = —1/4 and 
the model works for all time < t < oo provided the condition 2|2u;+3|0i > pi/a± 
is satisfied. The value of u is related to the other constants through the relation 



2cj + 3 



Pi 



a?0i 



± 



\ 



2af0! 



3pi 



6aiVi0i 



model works only for t > [ 2a a^ +3 | ] 3 - 4 



For other values of a the model does not work for the whole range of time 
0<t<oo. Ifa>4/3, when the rate of acceleration is faster than q = —1/4, the 

For 1 < a < 4/3, the model is valid 

only up to a time t = [ 2a a^2^+3| ] 3a " 4 during which the universe expands with an 
acceleration with a rate less than q = —1/4. 

For q = —1/4, where the model works for the entire time span, the present 
age of the universe can be calculated from ( El) as 



2V2 



\2u + 3) - 



PO 



ViA 4 ' 



1/2 



3(3a; + 4) 



1 



(18) 
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In the large u limit, this equation reduces to 



- (19) 

where u is obviously a negative quantity. Choosing V\ = w^ -1 and Ho ~ 65 
km-s _1 /Mpc, the present age of the Universe turns out to be approximately 15 
Gyr. 

Likewise, the present rate of variation of G is \G/G\q = \<p/<p\o = §#o < 10~ 10 per 
year. This is quite compatible with the observational data (see J7j and references 
therein) . 

All these possibilities are for (3 = —2. For other values of (3 also, one can obtain 
consistent solutions. Such as if (3 = — 1, the equation ( [14]) yields two possible 
solutions for a, namely 1/2 and —u/2. We disregard a = 1/2 as we are interested 
in non-decelerating models where a > 1. For uo = —2, we have a = 1, i.e. an 
uniformly expanding universe with q = 0. In this case, should be greater than 
pi/af. The equation system can be easily solved to get ip oc t -1 / 2 and V oc tp 6 . 
In this uniformly expanding scenario, the model can work for the whole range 
of < t < oo. If C(j is further negative, i.e., u < —2, we get faster rates of 
acceleration, but the model does not work for the whole range of time. 

Case 2 : k ^ 

In this case the condition ( |TJ]) remains in place and a further condition from 
IT]) and ( pD is 



(3 = -2a. (20) 
It can be seen from the field equations ( 0) and ( R]) that one can obtain 



consistent solutions only for a = 1. So equation ( |20[) immediately yields (3 = —2. 
With this value of (3, the condition ( [TJD is automatically satisfied for all values 
of uo. In this case also the model works for a limited period of time, < t < t±, 
where 



h = 20! 



\-(2cu + 3) 



^. (21) 
Pi 



(Xi 

For an open universe, i.e. for k = —1, (2uj + 3) has to be negative and 
\2lu + 3\ > l/a\. 

For a closed universe (k = 1), the model works even for a positive (2uj + 3) pro- 
vided (2uj + 3)af < 1. If however, (2u + 3) is negative, the model holds good 
without any such relation between a\ and u>. 

Thus, unlike most of proposed models with a non-positive definite deceleration 
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parameter, a quintessence field in Brans-Dicke theory works for a spatially non- 
flat Robertson Walker spacetime as well. It is true that for k 7^ cases the only 
consistent solutions have q = 0, but it is anyway non-decelerating and thus the 
expansion rate is faster than t 2 / 3 , which might sufficiently explain the recent ob- 
servations on the distant supernovae ||. 

It is worthy of note that the field equations can be integrated to produce consis- 
tent solutions for other forms of normal matter. Such as for a radiation dominated 
universe (p m = p m /3), a simple choice of V as V = VbV> 6 and a negative uj will 
lead to a decelerated expansion (q > 0) for the universe with a oc t* and oc t" 1 . 
Thus the model can be interpolated back to earlier epoch to yield a decelerated 
universe which is required in order to explain processes like nucleosynthesis. 



3 A possible solution to the flatness problem 

One important aspect of this model is that potentially it can solve the flatness 
problem as well. To see this we effect a conformal transformation as 

9n» = <P9iw, (22) 

which enables us to identify the energy contributions from different components 
of matter very clearly. Equation ( 0), in this new version, looks like 

3(a 2 + k) 



a 2 



Pm + P<f> + P^i (23) 



where a bar indicates quantities in the transformed version; p m , p^, p^ are the 
contributions to the energy density from the normal matter, the BD scalar field 
and the quintessence scalar field respectively. The quantity p^ is actually given 
by 



(2u + 3) 



(24) 



The quantities describing the matter distribution are transformed as pi = 
<fi~ 2 Pi and pi = <f)~ 2 pi. Now we define the dimensionless density parameter Q as 

^ = 3ll 2 = ^ m ^ 
where the total density p = p m + p^ + p^p, and the individual density paramaters 
Qi are defined accordingly. Using the equation ( ^3]) and the equation for the 
conservation for the total energy, 

p- + 35(p + p) = 0, (26) 
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we can write down the evolution equation for the density parameter as 

fl(fl-l)(37-2)# = 0. (27) 

The net barotropic index 7 is defined as 

7O = 7 m f2 m + 7^0^ + 7^^^. (28) 

The individual ^[s are defined by the equation pi — (73 — l)pi- The ratioes 
Pi/ Pi remain the same in both frames and thus 7j's do not vary. For our choice 
of matter, p m = and p^ = p^, and thus j m = 1 and 7^ = 2. The third 
index 7^, however, is not a constant and evolves with time via the equation 

li> = (jPi> + PtfVW = ^7(l^ 2 _+ V). 

Equation ( ^7|) indicates that Q — 1 is indeed a solution. The stability of this 
solution demands that (dti/dCl) H should be negative at Cl = 1. Equation ( p7]) 
shows that for an expanding universe (^ > 0) this is possible only if 7 < 2/3. 
From equation ( |28|) it can be shown that the relevant condition for 7 < 2/3 is 

ti m + 4fy < (2 - 37^, (29) 
which can be achieved by a suitable adjustment of the parameters. 

It is wellknown that geodesic equations are not valid in this conformally trans- 
formed version || and hence different quantities are not dependable regarding the 
content of their physical meaning. But it must be emphasized that the character 
of k remains unaltered, and thus if Qk = k/a 2 is zero in one frame, it must be so 
in the other as well. 

4 Concluding remarks 

A quintessence scalar field in Brans-Dicke theory is shown to give rise to an ac- 
celerated expansion for the present universe. Bertolami and Martins H] modified 
Brans-Dicke theory by introducing a potential function V = V{4>) where is 
the Brans-Dicke scalar field. As Brans-Dicke theory by itself is, in a sense, self- 
interacting (the kinetic term in the action contains 0), we do not include such 
a potential. Rather a non-gravitational field ip with a potential V = V(ip) is 
included. For a spatially flat universe, the model yields various non-decelerating 
solutions including a uniformly expanding solution (q = 0). For a simple choice 
of the potential (V oc ip A ), the model gives the solution of p. In this last case, 
the model can work for all time < t < 00. In most of the other accelerating 
solutions, the model is seen to work for only a restricted period of time. 

An important merit of this ansatz is that it provides solutions for a non flat 
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(k ^ 0) Robertson Walker metric as well. Although these solutions are not ac- 
celerating, they are not decelerating either (q = 0). So along with providing a 
non-decelerating solution, it can potentially solve the flatness problem too. In 
fact it has been shown that = 1 could be a stable solution in this model. 

This acceleration of the universe is achieved, in general, by a negative to. It 
has been claimed that the value of to should be large (> 500) if Brans-Dicke the- 
ory has to be consistent with the astronomical observations fllP] . But this value 
actually refers to the magnitude of to which can still be very large in this model 
except in some cases, such as for (k = 0,q = 0), where u = —2. Furthermore, 
reconciliation with Kaluza-Klein theory or low-energy string theory favours a 
negative value of u (see e.g. |Tl|). 

Acknowledgements 

This work has been partially supported by Spanish Ministry of Education under 
Grant PB94-0781. One of us (NB) is grateful to the "Direccio General de la 
Recerca" of the Catalonian Government for financial support under grant PIV99. 
The authors are grateful to Winfried Zimdahl for critically reading an earlier draft 
of this paper and useful suggestions. 

References 

[1] Perlmutter S et al. 1999 Astrophys. J. 517 565; Riess AG et al. 1998 Astron. 
J. 116 1009; Garnavich PM et al. 1998 Astrophys. J 509 74 

[2] Caldwell RR, Dave R and Steinhardt PJ 1998 Phys. Rev. Lett. 80 1582; 
Turner MS and White M 1997 Phys. Rev. D 56 R4439 

[3] Zlatev I, Wang L and Steinhardt PJ 1999 Phys. Rev. Lett. 82 896; P.J. 
Steinhardt, L. Wang and I. Zlatev 1999 Phys. Rev. D 59, 123504 

[4] Bucher M and Spergel D 1999 Phys. Rev. D 60 043505; Bucher M 1999 
[|astro-ph/ 99080471 

[5] Chimento LP, Jakubi AS and Pavon D 2000 Phys Rev. D (in the press) 
[|astro-ph/000507q| 

[6] Bertolami O and Martins PJ 2000 Phys. Rev. D 61 064007 

[7] Weinberg S 1972 Gravitation and Cosmology (New York: Wiley) 

[8] Riess AG 2000 "The case for an accelerating universe from supernovae" 
[|astro-ph/0005229| . 



8 



[9] Dicke RH 1962 Phys. Rev. 125 2163 

[10] Will CM 1993 Theory and Experiment in Gravitational Physics (Cambridge: 
Cambridge University Press, 3rd edition) 

[11] Kolitch SJ and Eardley DM 1995 Ann. Phys. (N.Y.) 241 128 



9 



